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THE SHRINKING FIGURE EIGHT AND OTHER SOLITONS FOR THE 

CURVE DIFFUSION FLOW 

MAUREEN EDWARDS, ALEXANDER GERHARDT-BOURKE, JAMES MCCOY*, GLEN WHEELER, 

AND VALENTINA-MIRA WHEELER 


Abstract. In this article we investigate the dynamics of special solutions to the surface 
diffusion flow of idealised ribbons. This equation reduces to studying the curve diffusion flow 
for the profile curve of the ribbon. We provide: (1) a complete classification of stationary 
solutions; (2) qualitative results on shrinkers, translators, and rotators; and (3) an explicit 
parametrisation of a shrinking figure eight curve. 


1. Introduction 


The curve diffusion flow is a particular evolution law governing the changing shape of a curve in 
time. It is the one-intrinsic-dimension analogue of the surface diffusion flow for evolving surfaces. 
These flows are gradient flows of a certain energy functional in a certain function space, but not 
the more common L 2 -gradient flow such as the curve shortening flow (corresponding to decreasing 
the length functional for a curve), the mean curvature flow (decreasing the area functional for a 
surface), or the Willmore flow (decreasing the elastic energy for a curve or the Willmore energy 
for a surface). 

Surface diffusion flow was first described by W. W. Mullins in 1956 as a model for the develop¬ 
ment of grooving at grain boundaries of heated crystal structures Mu . The flow is also related 


to the Cahn-Hilliard equation, which describes phase separation and coarsening in the quench¬ 


ing process of binary alloys CENC . For further analysis and historical remarks on the surface 
diffusion flow and its constrained variations we refer the interested reader to MWW Wh2|IWhl 


The surface diffusion flow is a fourth-order nonlinear parabolic system of equations that, at 
each time t, moves points p on the solution surface / perpendicularly to the surface with speed 
(A H)(p). Here A is the Laplace-Beltrami operator and H is the mean curvature, which is given 
by H = (A/, v). The inner product (•, •) is the standard dot product in R 3 . 

Surface diffusion flow is the geometric analogue of the classical parabolic clamped plate equa¬ 
tion, which expresses the dynamics of a flat plate realised as the graph of a function over a 
domain in R 2 . The curve diffusion flow is similarly the geometric analogue of the classical par¬ 
abolic Euler-Bernoulli model for a clamped rod, modelled as the graph of a function over an 
interval in R. When the domain is no longer flat, as is the case for a clamped arc segment of 
a circle in the plane or cylindrical ring in space, the geometry of the problem must enter into 


the model ESK SE . The surface or curve diffusion flow can capture the notion of stability for 
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Figure 1. Pieces of two cylindrical ribbons, immersed via f(u,z) = 
(x(u),y(u),z) where (u,d) £§xl. The curve 7 (u) = (x(u),y(u)) is called 
the profile curve. The ribbon on the left is embedded, while the ribbon on the 
right is immersed. The profile curve of the ribbon on the right is the self-similar 
figure eight curve we introduce in Section 6. 


an elastic plate or rod under stress respectively. Indeed, if the flow returns all small perturba¬ 
tions of a given initial configuration to that initial configuration without exception, then this 
configuration is called stable |Nj. 

If the evolving surfaces are cylindrical, as would be the case for idealised closed ribbons with 
their width extended to infinity, then this remains so under the flow. In this case the flow 
equation reduces to a one-dimensional parabolic equation for the profile curve 7. For the surface 
diffusion flow, this equation is exactly the curve diffusion flow. This reduction is also performed 


SE 


Let us briefly remark on the above procedure with respect to the particular case of a Mobius 
band. Extending the width of the Mobius band to infinity must result in self-intersections, and 
so, the profile curve does not project uniquely to a plane curve in R 2 . This is an essential 
characteristic of the Mobius band. One way in which we could overcome this difficulty is instead 
to immerse the Mobous band in R 4 , then extend the width to infinity, and project to R 3 . This 
would produce a profile curve in R 3 , the analysis of the dynamics of which are much more 
complicated than for the planar curve case. 

In this paper we concentrate on the dynamics of solutions to the curve diffusion flow which 
evolve either by scaling, translation or rotation only. These solutions are known as solitons. 
Solitons are of interest both from a theoretical and practical viewpoint due to their relationship 
with singularities: given a family of evolving curves, if the family exists for at most a finite amount 
of time, then it must develop a singularity at the final time. Upon rescaling the flow, effectively 
‘zooming in on the singularity’, we expect to construct a family of curves which approximate one 
of the solitons for the flow. 

The existence of finite-time singularities is typically a difficult question, especially in the case 
of fourth and higher-order curvature flows. For the curve diffusion flow, this was settled in 
Polden iP . The analysis there was later extended by Escher and Ito 
configurations. In particular, they prove the following theorem. 


El to more complex initial 


Theorem 1 (Escher-Ito-Polden). Let 7 0 be a plane curve with winding number zero. Let T 
denote the maximal time of existence for the curve diffusion flow 7 1 := 7(71) with initial data 
7o = 7(‘,0). Then 
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(i) The estimate T*("f 0 ) := > T holds, where L( 7 0 ) denotes the length of 70 ; and 

(ii) If T = then the solution shrinks to a point as t approaches T. 

This theorem raises two interesting questions: 

(1) Is there a solution with winding number zero whose maximal existence time satisfies 
t = t *( 7 )? 

(2) Is there a solution with winding number zero that shrinks to a point? (That is, no other 
kind of singularity occurs.) 

The first question investigates the sharpness of the theorem. Its answer is no, since in the proof 
of Theorem |T| Fenchel’s inequality is used, and for plane curves with winding number zero, it is 
strict. Furthermore, the upper bound from Theorem 1 may be improved by a factor of yy by 
using the following argument]]] 

Let 7 1 be a family of curves evolving by curve diffusion flow. Differentiating the length yields 

( 1 ) X'( 7 f) = - f \d s n\ 2 dl, 

J-f 

where k(s) is the curvature of 7 t and 7 (f) is parametrised by arc length. We use d s and d£ to 
denote the arc length derivative and element respectively. (We refer the reader to Section 2 for 
more detail on the notation.) For curves with winding number zero, the curvature n has zero 

average, and furthermore ^ J] , M ddj <f^L(7i)/ 7 l «| 2 d£ (this holds for any function with zero 
average), so, using the standard sharp Poincare inequality we compute 

( l w dt ) -B L l M 2dl -f. L ^l =3s l la “ |2<K - 

In the above we have omitted the argument 7 * from L. Since ||k||i := f |/c| d£ > 2-7T (by a simple 
variant of Fenchel’s theorem), we combine this with |l]) above to conclude 

(L 4 )' = 4L 3 L' = -4||<9 s k|| 2 L 3 < -4(48||re||?) < -3 • 4V 
=> L 4 ("/ t ) < -3 • 4 4 7 r 2 t + L 4 ( 7 o) , 

and consequently we obtain the estimate T < T("/ 0 ) := L 4 ( 7 0 )/ 7687 r 2 , which is better than the 
estimate given by Theorem 1. 

In the above proof, there are three key estimates used: Fenchel’s inequality, the Poincare 
inequality, and a sharp Holder-type inequality for functions with zero average. It is not the 
case that the same curve satisfies equality for all three of these and therefore it appears clear 
that the improved bound T < T( 70 ) is again not optimal. The broader question therefore on 
characterising the optimal upper bound on the maximal time of existence in terms of initial 
length for curves with winding number zero remains open: 

(1’) Is there a T depending only on initial length L( j 0 ) such that: 

(a) All curve diffusion flows with winding number zero have maximal time bounded by 
T; and 

(b) Does there exist at least one curve diffusion flow with maximal time equal to T, or, 
does there exist a sequence of solutions 7 j with initial data 7 ] 0 * and winding number 
zero whose maximal existence time X) satisfies |X( 7 f°' ) ) — T)| — > 0? 

In this paper we shall not answer this question. We do however provide an answer to the 
second question (2) above. The second question is interesting from a classification of finite¬ 
time singularities perspective. The key issue there is to demonstrate that solutions may become 


4 We would like to thank one of the anonymous referees for pointing out this improved bound. 
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singular by simply becoming small - a rescaling of such solutions will exist for all time and be 
asymptotic to a smooth model curve. A finite-time singularitiy of this type is said to be of Type 
I, in analogy with existing work on the mean curvature flow Hu . Other singularities are called 
Type II. Type I singularities are expected to be simpler in structure and satisfy natural stability 
criteria when compared with Type II singularities. 

We answer the second question in the affirmative by way of an explicit example: a self-similarly 
shrinking immersed figure eight. This figure eight is in fact the lemniscate of Bernoulli^] Since 
one requires zero area to have the possibility that the solution shrinks to a point (see Lemma 
[5]), and the maximal existence time is finite (see Theorem 1), a figure eight with one leaf smaller 
than the other will certainly become singular with non-zero area, ruling out the possibility that 
it is asymptotic to a point. 

Let 7 be the shrinking figure eight described in Section 6 . Upon scaling 7 to satisfy L(y) = 1, 
we compute the extinction time of the flow to be 


3 • 2 11 ■ DC(-l ) 4 ’ 

where 3C(m) is the complete elliptic integral of the first kind with parameter m = k 2 . In this 
particular case the maximal time estimate T* overshoots T by a factor of roughly 2.9, and T 
overshoots T by a factor of roughly 2.4: 

rp* rp 

— = 2.9115257845 ... , - = 2.3946339747.... 

In a sense, the self-similar figure eight represents a “best-case” scenario. The generic case is that 
an arbitrary smooth curve with winding number zero will be driven to a curvature singularity 
and not shrink to a point. Furthermore, such a singularity involves blowup of the curvature 
and its derivatives, so in particular the speed of the flow is greater in this circumstance. These 
heuristics lead us to make the following conjecture. 

Conjecture. The maximal existence time of a curve diffusion flow emanating from a curve with 
winding number zero is bounded by the existence time of a self-similar figure eight with the same 
initial length. 


Of course, verifying the conjecture has uniqueness of the self-similar figure eight as a corollary. 

Another popular fourth-order surface evolution equation is the Willmore flow, whose one¬ 
dimensional counterpart, the elastic flow, also models cylindrical elastic ribbons. The key dif¬ 
ference between the two flows is that the elastic flow minimises the integral of the square of the 
curvature of the curve, while the curve diffusion flow minimises length, keeps enclosed volume 
fixed, and keeps the quantity 


«-W, 


\d s K\ 2 dt ) dt 


uniformly bounded, with a bound dependent only on the initial data. A proof of these statements 
can be found in Wh3 . We remark on the similarities and differences between this flow and the 


curve diffusion flow throughout this article. 

The structure and specific contributions of this article are as follows. We begin Section 2 
by setting up notation, defining key geometric quantities and giving the fourth-order nonlinear 
parabolic system of partial differential equations for the curve diffusion flow. We then briefly 
review relevant previous work and the question of local existence before moving on to a clas¬ 
sification of stationary solutions - those that remain constant under the flow. We provide a 
five-parameter family of parametrisations which completely classifies these solutions. In Section 


2 The authors would like to thank Hojoo Lee for pointing this out. 
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3 we record some elementary properties of the flow that are relevant to our investigations. In 
Section 4 we focus on curves that evolve self-similarly, deriving the corresponding fourth-order 
ordinary differential equation satisfied by such curves. We prove several properties of self-similar 
solutions. Section 5 is devoted to curves that evolve by translation only. We prove there that 
the only closed translators are circles, and that the only open translators (open immersions of R) 
satisfying either of two conditions at infinity are straight lines. In Section 6 we provide an explicit 
parametrisation of a figure eight immersed curve and verify that this curve does indeed evolve 
self-similarly under the curve diffusion flow, contracting to a point in finite time. In Section 7 we 
consider the case of rotating solutions and prove some non-existence results analogous to those 
for translators in Section 5. 
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2. The curve diffusion flow equation 

We will parametrise closed evolving curves over the unit circle S, and open evolving curves 
over the real line R. Let 7 denote the evolving curve, with components 

7 (u) = {x (u),y(u)). 

The tangent direction to the curve is given by 

lu (u) = ( x u ( u),y u (u )), 
and therefore a choice of unit normal is 

v (X) = j— ^77 (y u («), -x u (u )), 

\lu (u)| 

where |j denotes the length of a vector. We will always assume our curves are smooth and 
regular, that is, | 7 „| > 0 for all u. Recall that, in a general parametrisation, the arc length of a 
curve 7, beginning at point 7(110), is given by 

(2) s(u) = f | 7 „ (u) | du. 

Ju 0 

The arc length induces an arc length element d£ = | 7 „| du. The curvature vector and curvature 
of a plane curve is given by 

k = 7 ss and k = {is, k) , 

where we use (•, •) to denote the standard inner product between vectors in the plane. For closed 
embedded curves, we will always choose an inward-pointing unit normal, so that the curvature 
of a circle is positive. Otherwise, we select either normal - the curve diffusion flow is invariant 
under changes in orientation. 

We say = j(-,t) evolves under the curve diffusion flow if it satisfies the system of fourth- 
order partial differential equations 

0 'y 

(3) -^(u,t) =-K as {u,t)v(u,t) for (w,t) £§x (0, T) 
with initial condition 


7 (u, 0 ) = 70 (u), 



6 


M. EDWARDS, A. GEHARDT-BOURKE, J. MCCOY, G. WHEELER, AND V.-M. WHEELER 


for some given initial curve immersion 7 . Note that an immersed curve may have intersections, 
however, it has a well-defined tangent vector for each u. We label the maximal existence time of 
a solution to ([3]) as T. 

We are interested in the shape of the evolving curves, that is, the image 7 ( 8 , t) for closed 
solutions, or 7 (R,f) for open solutions, as a set in M 2 , and not in the orbit of a specific point. 
There is a degeneracy in the flow equation ([ 3 ]): reparametrisation and tangential movement of 
points in the direction of 7 S do not affect the shape of the image 7 ( 8 , t). The degeneracy consists 
of every diffeomorphism <p : 8 —>■ 8 . This implies that, in terms of the images 7 ( 8 , f), the flow 
0 is equivalent to 

(4) (u,t) =-K ss (u,t) for (u,t) e 8 x (0,T) 


The impact of this observation on our analysis here is that when considering solutions with a 
particular form, we must also take into account a possible diffeomorphism p acting on 8 . This 
typically forces us to take normal projections in order to remove the influence of p. 

Note that 


9s2(7ss,i/) hd ( k,l ( (l'Yj(hd) 


1 


1 


l7ul \l7ul V \ I 7 J MTu 

so that the leading order term of — k ss = v ) is precisely 



1 


\lu 


(lu 


■') • 


The normal direction is ‘positive’ for the flow ([ 3 ]), and so this shows that ([3| is a quasilinear 
parabolic systenj^J 

Current understanding of the curve diffusion flow is quite limited, apart from short time 
existence. There are several methods for obtaining short time existence of solutions to ([3]) given 


suitable initial data; we refer the reader to the discussion of these in Wh3 and the references 


contained therein. Very little else is known about the behaviour of the flow. Giga and Ito (Gil 


provided an example of a simple, closed, embedded plane curve that develops a self-intersection 
in finite time under the flow @. They also gave an example in [GI2| of a closed, strictly convex 


plane curve that becomes non-convex in finite time. Elliot and Maier-Paape showed that curves 
that are initially graphical may evolve under ([3| to become non-graphical in finite time EMP 


On the other hand, the fourth author recently showed Wh3 that under the curve diffusion flow, 


solution curves that are initially close to a circle, in the sense of the normalised integral of the 
square of the oscillation of curvature, exist for all time and converge exponentially fast to a 
(possibly different) circle. Further, those admissible curves which are initially embedded remain 
embedded. A fundamental obstacle in the analysis of such flows and of ([ 3 ]) is that, as a system 
of higher order partial differential equations, arguments based upon maximum principles are in 
general not available. 

Let us now move on to the classification of stationary curve diffusion flows. We call the family 
of curves 7 (-,t) a closed solution if 7 satisfies 0 and 7 (u) = 7(it + P) for some P > 0. The 
smallest such P is the period of 7. 

A solution is called stationary if it does not change in time under Q, that is, if K ss (u,t) = 0. 

Lemma 2 . Suppose 7 : 8 —> R 2 is a closed stationary solution to the curve diffusion flow. Then 
7 is a round circle. 


:i Tt may be helpful for the reader to recall the biharmonic heat flow ( 0/ + A 2 )v = 0 for a function v. 











SOLITONS FOR THE CURVE DIFFUSION FLOW 


7 


Proof: Integrating by parts and using k ss = 0 we have 

0 = — I k k ss (U = / Kj d£ 

J 'y J / y 

and so k s = 0, that is, n is constant. Therefore 7 is a smooth closed curve with constant 
curvature, and by the classification of curves in IR 2 , any curve with constant curvature in the 
plane is a straight line or round circle. Since 7 is closed, 7 can not be a straight line, and so it 
must be a round circle. □ 

The situation is more complicated for open solutions. We define an open curve to be a regular 
immersion 7 : R —> IR 2 of 1 R for which there does not exist a pair of intervals /, J such that 

( 5 ) In J = 0 , |/| = 00 , |J|<oo, and 'y(I) = j(J). 

This last condition captures the notion of not allowing 7 to ‘close’ while still allowing infinitely 
many self-intersections, including possibly intersecting on an arbitrarily large open set. 

Open curves may concentrate around a particular point, such as a logarithmic spiral. This 
concentration behaviour is encapsulated by the notion of properness of the curve 7 : ]R —> K 2 . 
The curve 7 is proper if for any compact set CcR 2 the inverse image 

7 _ 1 (C) := {«el: 7 («) G C} 

is a compact set. 

Lemma 3 . Suppose 7 : 1 R —» ]R 2 is an open stationary solution to the curve diffusion flow. Then 
7 is either a straight line, or a similarity transformation of the standard Cornu spiral. Up to 
translation and rotation such curves satisfy 

k(s) = fes + jfei 

for a pair (k\, G 1 R 2 . 

Proof: If k ss = 0 , then there is a &2 G R. such that 

k s (s) = k 2 . 

The curvature n is therefore linear in s; that is, there exists a k\ G K. such that 

k(s) = lt2S + jfei. 

If /c2 7^ 0 , then after reparametrising by t = s + we have 

K(t) = k 2 t ■ 

Furthermore, by scaling the curve 7 by a factor of \fk2l we have 

k{t) = ^ r — K,(t) = \[kzKt = 7 it. 
v ^2 

Therefore either k 2 = 0, in which case we have a straight line, or k2 ^ 0 and we have a similarity 
transformation of the standard Cornu spiral. □ 

Remark: Combining the above two lemmata, we see that up to similarity transformations the 
only three stationary solutions are the circle, straight line, and Cornu spiral. 

Corollary 4 . Suppose 7 : 1 R —► K 2 is an open, stationary, properly immersed solution to the 
curve diffusion flow. Then 7 is a straight line. 
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Proof: We apply Lemma [ 3 ] above to obtain that 7 has curvature satisfying 

k(s) = k 2 s + k \ . 

If k 2 = 0 then k is a constant and so 7 is a straight line, as required. 

Suppose on the contrary that k 2 7^ 0 . We will show that this leads to a contradiction to 
the properness assumption. Consider the restriction 7 of 7 to the interval ( 0 , 00) if k 2 > 0 or 
to the interval (—00, 0 ) if k 2 < 0 . The curve 7 has strictly monotone positive curvature. We 
may therefore apply the classical Tait-Kneser theorem to obtain that the osculating circles of 7 
are pairwise disjoint and nested (see fGTTl and the references therein for a discussion of Tait’s 
original paper and some interesting extensions). Recall that the osculating circle at 7(so) is the 
best approximating circle to j(so). It is a standard round circle with radius 1 / k ( so ) and centre 
(7(so) + v( s o)/k(s 0 )). 

Let C(so) denote the osculating circle at 7(so) and D(s0) denote the disk with boundary C(so)- 
The theorem implies that 

®(s 0 )D (J C ( s)D (J 7(a). 

SG[so,Oo) sE[so,OC>) 

Clearly D(so) is compact. The above computation shows that 

7^ 1 (B(s 0 )) = [s 0 ,00), 

which is not compact, and so 7 is not proper. This is contrary to the hypotheses of this Corollary, 
and so k 2 = 0, as required. □ 


The two-parameter family of open curves given by Lemma [3] therefore consists of a particular 
family of improperly immersed curves and straight lines. The improperly immersed curves are 
Eulerian spirals, with the clothoid or Cornu spiral (k 2 = 7 r and k\ = 0) a standard example. The 
arc length parametrisation of the clothoid is given in terms of the Fresnel C and S integrals: 

7 (s) = ( C(s), S(s) ), where C(s) = J cos dt and S(s) = J sin dt. 

The curvature of the clothoid is k(s) = S7r. In order to produce the whole two-parameter family, 
let us define the modified Fresnel integrals C(s, c\,c 2 ) and S(s,ci,c 2 ) by 

C(s, Ci,c 2 ) = / cos (cit + c 2 t 2 ) dt and S'(s,ci,C 2 )= / sin (cit + c 2 t 2 ) dt. 

Jo Jo 

The curve 7 : R —► K 2 with arc length parametrisation 

7 (s) = ( C(s, ci, c 2 ), S(s, ci, c 2 )) 


then has curvature k(s) = 2sc2 + Ci. 

Note that this includes the case of circles (C2 = 0 ) and straight lines (ci = C2 = 0 ). The 
curvature k is invariant under rigid motions (rotations and translations), and therefore in order 
to account for the entire family of stationary solutions we must include all possible rotations and 
translations of the model members y(s) = (C(s,ci,c 2 ), S(s,c\,c 2 ) ). Therefore the family T of 
immersed open or closed stationary solutions to the curve diffusion flow consists precisely of the 
five parameter family 


T = 



cos 0 — sin 6 

sin 6 cos 9 


C(s, Ci, c 2 ) 
S ( s , C !, C 2 ) 


+ (V U V 2 ) for some (d, c 2 , Vi, V 2 ,9) e R 4 xS 


up to reparametrisation. Some members of T are depicted in Figure [ 2 ] 
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Figure 2. Three members of T: the left with k(s) = —2s + 3, the middle with 
k(s) = s , and the right with k(s) = 4s + 13. We choose { 6 , V\, V?) = (0,0, 0) for 
each. 


Remark (Elastic flow): We say 7 1 evolves under the elastic flow if 7 (-,f) = 7 t satisfies 


d'y 

dt 



- K ss (u,t) 



The elastic flow is the steepest descent L 2 -gradient flow of f k 2 dl. For global regularity results 

Stationary solutions are termed elasticae and 


and further historical references we refer to 
satisfy 


DKS 


Kss{u,t ) = -~K 3 {u,t) . 


A classification of solutions to this equation is not yet known in the sense discussed here, that 
is, in terms of explicit parametrisations. Indeed, there are no known closed solutions to this 
equation. Open solutions are yet to be classified. We conjecture that there are no open properly 
immersed solutions apart from the straight line. 


3. Elementary properties of the curve diffusion flow 


In this section we examine the behaviour of the arc length of 7 1 , and the area enclosed within 
7 t, under the flow <©• Because 7 1 may be an immersed curve, care is required in defining the 
enclosed area, or simply the change in enclosed area; we refer the reader to the discussion in 


MW for details (that article concerns evolving immersed surfaces and discusses the definition 


of the signed enclosed volume, but the issues are analogous). When the curve 7 * is an embedded 
curve, the calculations below are the one intrinsic dimension analogues of those in Mcl , for 


example. The full details for the specific case of the curve diffusion flow of immersed curves are 
given in 


Wh3 


Lemma 5. Under the flow ([3]), the area enclosed by the curve 7 1 remains constant. In particular, 
if the family of curves contract to a point, then the initial data 70 must have zero signed area. 


Proof: The signed area A (jt) enclosed by evolves under ^ according to 

^ A (7 1 ) = - j k ss d£ = 0 . 

That is, A (t) is constant under the flow. 

Since the area of a point is zero, and the enclosed area is constant under the flow, the area at 
any time along a family of curves that contract to a point is zero. In particular, the area of the 
initial data 70 must be zero. □ 


Lemma 6. Under the flow ([ 3 ]), the length of the evolving curve 7 1 is monotone nonincreasing. 
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Proof: Under the flow (|3|, the length of the evolving curve 7 evolves according to 


d_ 

dt 


L(t ) 


K ss K.d£ = — k 2 s dt < 0 , 

< J -y 


that is, the length is nonincreasing. 


□ 


The above lemmata make the curve diffusion flow interesting from an isoperimetric point of 
view. In particular, the isoperimetric ratio is monotonically decreasing: 

Corollary 7 . The signed isoperimetric ratio 


Ii.lt) = 


L\lt) 


47rA(7 t ) 

is decreasing in absolute value with value given by 

I("/ t ) = /(7o)e _ ^° i(-rt) S-,*' dt * r . 
Proof: Differentiate I and apply Lemmata [5] and [ 6 ] to obtain 

d T( , 2 lint) f 2 


Integrating yields the result. 


□ 


Remark (Elastic flow): Although the elastic flow reduces f k 2 dt with velocity 


dtl/ M ~ 2 


k ss + 2 K " 


dt, 


the behaviour of the length, area, and isoperimetric ratio, are for generic initial data unknown. 


4. Curves evolving homothetically 


A self-similar solution to a curvature flow equation is a solution whose image maintains the 
same shape as it evolves, that is, it changes in time only by scaling, translation and/or rotation. 
In this section we focus on the case where 7 evolves by scaling, and in the next section we study 
solutions which evolve by translation. Section 7 studies curves evolving by rotation. 

Closed, self-similar solution hypersurfaces of positive mean curvature for the mean curvature 
flow are known to be spheres Hu ; the third author generalised this result for a class of fully 
nonlinear curvature flows in Mc2 . In the corresponding case of the curve shortening flow, 


self-similar curves were classified in |AL|. 

Definition. We say 7 1 is a solution of ([3]) (or equivalently Q) evolving homothetically if and 
only if there is a differentiable function of time / and an initial curve 70 such that 

7 (u, t) = f (t) 70 (u) . 


This is equivalent to 

(d t j,n) = f (70, v) . 

In general, homothetic solutions may be expanding, shrinking, stationary, or exhibit more com¬ 
plex breathing behaviour. For the case of curve diffusion flow, Lemmata [5] and [ 6 ] show that any 
solution must decrease its length monotonically (and be asymptotic to a curve with constant 
curvature unless it exists only for a finite amount of time). We therefore term a homothetic 
solution 7 a shrinker or shrinking solution. We take / (0) = p where p £ R as initial condition. 
This freedom will be useful when we estimate the maximal existence time of the shrinking figure 
eight in Section 6 . 

The simplest example of a closed shrinking curve which is a solution to (|3j) is the circle. 










SOLITONS FOR THE CURVE DIFFUSION FLOW 


11 


Lemma 8. Let 7 be defined by 7 (it, t) = R (cos u, sin it) for some constant R > 0 . Then 7 is a 
solution to (| 3 |. 

Proof: We simply observe 

Q 7 Q 

(it, t) = —R (cos 11, sin it) = 0 
dt y dt y ' 

and since the curvature k of 7 is identically equal to 

k ss (it, t) = 0. 

Hence © is trivially satisfied. □ 


If a curve is evolving homothetically under ([3]), we may write down a corresponding ODE by 
‘separation of variables’ that it must satisfy. 

Lemma 9 . Let 7 : S —> R 2 be a curve. Then 7 is a self-similar solution to ([ 3 ]) iff there is a 
constant K such that 

(6) k ss (it) = K (7 (it), v (it)) 

for all 11 £ S. 


Proof: Suppose the curve 7 is such that 7 (it, t) = f (t) 70 (it, t) for some positive differentiable 
function / (t) with / ( 0 ) = p , that is, at the initial time t = 0 , 7(11,0) = P7o (it). Geometric 
quantities associated to 7 are related to those of 70 as follows. In the calculations here we often 
omit the arguments of 7, 70 and /. The tangent direction is 

7u = f dulo 


and the unit normal is 


v = v 0 . 

Above and in what follows we will include a zero subscript on geometric quantities associated 
with 7 q. The curvature of 7 is related to that of 70 by 

(7) K = f- 1 K 0 . 

In view of ©, the arc length of the curves 7 and 70 are related via 

s = fs 0 

and so differentiating Q twice, we find 

Has = f~ 3 d 2 SoSo K 0 . 

Using the above we find that (J 3 J) is satisfied if and only if 

/' (7o, u 0 ) = -f~ 3 d 2 QSo K 0 (u 0 , v 0 ) ■ 

Following a standard separation of variables argument (note that any zeros of (70, v 0 ) are isolated 
unless 7 is a straight line passing through the origin - this case is trivial and can be treated 
separately), there must be a constant K such that 

(/ 4 )' = - 4 ^ SoS ° K ° = 4 K . 

\lo,vo) 

Simplifying this equation finishes the proof. □ 


Remarks: 



12 


M. EDWARDS, A. GEHARDT-BOURKE, J. MCCOY, G. WHEELER, AND V.-M. WHEELER 


1. Solving the ordinary differential equation for / in the above proof, we have 

/ (t) = {/p^ + AKt. 

Assuming also that © is satisfied, we observe that 

• If K = 0, then f (t) = 1 and the curve 7 is static under ©, that is, a stationary 
solution that exists for all time t. 

• If K < 0, then the solution to © shrinks self-similarly to a point at time T = — . 

• If K > 0 then any self-similar solution would have to expand indefinitely. However, 
for closed curves this is not possible in view of Lemma [ 6 ] For open curves this 
is possible, and there could be open expanding solutions. These could be impor¬ 
tant models for how the curve diffusion flow smooths out isolated singularities, in 
the same way that the error function generates a self-similar solution to the heat 
equation with the Heaviside step function as initial data. 

2. Clearly, the shrinker condition © is invariant under translation and rotation; that is, 
applying these operations yields another solution to ([ 6 ]) with the same constant K. Un¬ 
der scaling by a factor p, one again obtains a solution to © , but with a new constant 
K = Kp 4 . One may express the ODE © in terms of a pair of equations for (x(u),y(u)) 
in the original (non-arclength) parametrisation, and apply a standard Lie symmetry 
analysis to the system. The result of this analysis is that there are no other possible 
transformations that one may apply in order to obtain more solutions. 


Remark (Elastic flow): One naturally expects expanders for the elastic flow, in contrast with 
the shrinkers of curve diffusion flow. The solution to the elastic flow with initial data equal to 
a circle expands under the flow for all time, and so the standard round circle is the canonical 
example of an expander. Indeed, Theorems 3.2 and 3.3 of DKS establish global existence and 


convergence to an elastica under the condition that the length of the evolving curves remain 
fixed. Without this constraint, it appears impossible to obtain convergence of the flow in any 
sense, since there is no way to control the length. There are no other known expanders for the 
elastic flow. 


( 8 ) 


5. Curves evolving by translation 
A family of curves 7 : S x [0, T) —> R 2 evolving purely by translation satisfies 

cly 


— (u,t) =7o(u)+Ut + 7 s (w,t)0(u,t) 


for some constant vector V and smooth diffeomorphism cf>. In this case 

%•') = ^ • 

and if © is to hold then 7 must satisfy 

(9) - k ss = ( V , v). 

We call the solution 7 a translator. If V = 0, then the solution 7 is stationary, and we call the 
translator trivial. 

It is common that for a given curvature flow there are no closed non-trivial translators. Nev¬ 
ertheless one typically expects open immersed non-trivial translating solutions. A family of open 
immersed curves 7 : R x [0, T) —» R 2 evolving purely by translation satisfies © and © as in 


the closed setting. Such curves usually arise as blowup limits of a Type 2 singularity [GH 
Singularities of this type have engendered much interest in the literature. 
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In this section we investigate translators for the curve diffusion flow. We prove by a simple 
argument that there are no smooth closed non-trivial translators. Our proof implies a stronger 
statement about open curves, which is sharp, by the earlier example of the clothoid. 

Proposition 10. Let 7 : S —> M 2 be a smooth closed translator. Then 7 is trivial; that is, y(S) 
is a standard round circle. 

Proof: Since 7 is closed, integration by parts gives 


( 10 ) 


— / nn ss di = / n 2 „di. 


Also, noting that V is a constant vector and that the curvature vector ft satisfies nv = R = 7 ss , 
we have 


(ID 


[ (v,z) dt= (d s ({v, ls ))dt = 0 . 

J 'y J y 


Multiplying Q by k and integrating gives 

— / kk ss di = (V, , ft) dl , 

J'y J -y 

which when combined with (101 and 0 above yields 

j n 2 di = 0 . 

J 'y 

Since 7 is smooth, this implies that k s e 0 on 7 , that is, that the curvature is constant. Therefore 
ft ss = 0, and so V = 0; that is, the translator is trivial. By Lemma 2, it is a round circle. □ 

Now let us investigate open curves. 

Proposition 11. Suppose 7 : R —> M 2 is a smooth open translator satisfying 


( 12 ) 


lim ( K(pi)K s {pi ) + (V,7 s (pi)) ) - (ftfe)ft s (<?i) + (P,7 sidi)) ) 

1 —^00 L / 


= 0 


for arbitrary sequences {p^} and {qi} asymptotic to —00 and +00 respectively. Then 
straight line. 

Proof: Let us follow the idea of Proposition 10 We integrate ([9]) against ft to find 

r<H r<n 

/ n 2 „d£ = — / ftft ss di + ftft s ' 


is a 


(V,, ft) di + ftft. 


I Qi 
I Pi 


d s ( ( V, 7 S )) di + ftft. 


I Qi 
I Pi 


{V, 7s) + ««s 


J Pi 


Taking i —> 00 yields the result. Note that circles are excluded by the open hypothesis. 


□ 


Remark: The result of Proposition 11 is false without condition ( 12 ). To see this, consider the 
clothoid: a smooth open curve with ft(s) = s where s is the arc length parameter. Clearly this 
curve satisfies n ss = 0 and so is a stationary translator with V = 0 in ([£]). The clothoid however 
does not satisfy condition ( 12 ), since k(s)k s (s) = s is an odd function, and so the limit does not 
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exist. 


Although the above remark shows in some sense the sharpness of Proposition 0 it leaves 
open the possibility of replacing the decay condition ( 12 ) with something different. 

The following proposition explores this and converts (12) from a pointwise decay condition on 
the curvature to a natural growth condition on n in L 2 and P in L 1 . 

Proposition 12. Suppose 7 : R. —> R 2 is a smooth open proper translator satisfying 


(13) 


lim 


1 


k 2 d£ + 


p_> o° h~ x {B p { 0)) 


-f 

P d 7 -i(S 2 p ( 0 )) 


(V,7.)| dl.) =0, 


where B p { 0) denotes the closed ball of radius p centred at the origin. Then 7 (R) is a straight 
line. 

Proof: We use again a similar idea, however this time we localise the estimate in the plane using 
a cutoff function ?y p = fj p o 7 where fj p : R 2 —> [ 0 , 1 ] satisfies 

Vp(x) = 1 on Be{ 0), 

Vp( x ) = 0 on B p { 0), 
rj is of class C°° , supp ry is compact, and 

there exist constants Ck such that d^rj < . 

pk 


The existence of such a function is straightforward. A constructive proof can be found in Wi 
Let us denote by r )~ l {S) the inverse image of a set 5 C R 2 . First, note that since 7 is proper, 
the integral 

f re 2 ?y 4 dl 

d 7 -i(B p (0)) 

is well-defined, irrespective of p. We compute 


k 2 ry 4 d£ = — 


/ 7 - 1 (b p ( o)) 


/ 7 - 1 (b p ( 0 )) 

/ 

/ 7 - 1 (B p ( 0)) 


kk ss ry 4 dt — 4 


d 7 -i(B p (0)) 
= -4 


kk s rj s ?y 3 d£ 


(P, k) ry 4 dl — 4 f kk s ry s rf d£ 

i 7 -MB p (0)) 

d s ((P, 7 S )) ?y 4 dl - 4 f kk s ?y s 77 s d£ 

d 7 -i(B p (0)) 

/ {V^s) VsV 3 d£ - A kk s 7j s rf dl 

d 7 -i(B„(0)) J'y-UBJO)) 


< £ 


/ 7 - 1 (B p (0)) 

k 2 ?y 4 d£ + 4 


J 7 -i(B p (0)) 

4 f 

+ 


/ 7 -i(S p (0)) 


7“ 1 (Bp(0)) 

(Pj 7s) I |?y s |7y 3 ^ 


k 2 ry 2 ry 2 d£ , 


£ -V*(B p (0)) 

for any e > 0. Choosing e = | and absorbing the first term on the left, we find 


n 2 dl- 


/ n 2 s ifd£<^ 

J'l-HBpm P d 7 -i(B 2p (0)) P d t _ 1 (B 2p (0)) 

Taking p —>■ 00 in this estimate yields k s = 0 and hence the result. 


(V,y>)\dt. 


□ 
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Remark: It is an open problem to determine if either the properness condition or growth 

One advantage of the integral growth condition 


condition ( 13 ) are necessary in Proposition 12 


over the pointwise decay condition is that the rigidity statement holds for a much wider class of 
weak solutions, which may not even have continuous tangent vector. We will not explore notions 
of weak solution here. 


6. The shrinking figure eight 
Proposition 13 . Let 7 : § —>• R 2 be defined by 


( 14 ) 


7 (u) = 


1 


1 + sin u) 


Then 7 evolves self-similarly under <§• 


cos u , - sin 2 u 


Remark: The curve defined by (| 14 j) is the lemniscate of Bernoulli. 
Proof: The tangent direction to 7 is given by 

1 


In (u) = 


(2 — cos 2 u)‘ 


(— (2 + cos 2 u) sin u, — 2 + 3 cos 2 zt) 


and the unit normal is 


v(u) = 


3 sin 2 zt - 1 - sin u (2 + cos 2 u) 


Therefore, we have 
( 15 ) 


(l+sin 2 it) 2 (l + sin 2 zt) 2 
cos u (sin 2 u — l) 

(7 (it), v (zt)) = 


(l + sin 2 zt) 2 

By further straightforward calculations, the curvature of 7 is given by 

, . 3 cos (zt) 

« W =--—r- 

(l + sin 2 zt) 2 

We calculate using ([2]) with the chain rule that the derivatives or curvature with respect to arc 
length are 

. —6sinu — 6cos 3 u 

K s (u) = -r~2— an d K ss (w) 


I + sin u 


(1 + sin 2 zt)' 


Comparing with ( 15 ), we see that 

6 (7 (it), v (zt)) = k ss (zt), 
so 7 given by ( 14 ) satisfies © with K = —6. 


□ 


Remarks: 

1 . Since K < 0 above, in view of the Remark at the end of Section 4 , the figure eight given 
by ( 14 ) shrinks to a point at time T = ^. This is in fact consistent with Lemma pi 
since the symmetry of the figure eight implies that the signed area is equal to zero and 
remains constant under ([ 3 |. More generally, in view of Lemma [ 5 ] any curve that shrinks 
self-similarly under (| 3 j) must have zero signed area. 
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Figure 3. The shrinking figure eight (middle), along with plots of its curvature 
(left) and second arc length derivative of curvature (right). The horizontal axis 
represents the u variable and the vertical axis represents k(u) and k ss (u) in the 
left and right figures respectively. 


2. The length of the parametrisation of the figure eight above is 

L{ 7) = 43C(-1). 

Therefore, choosing p = 4 g C (_ 1 ) allows us to obtain the maximal existence time of the 
figure eight with unit length as claimed in the introduction. 


7. Curves evolving by rotation 

A family of curves 7 : S x [0, T) —> R 2 evolving purely by rotation satisfies 
(16) - K sa = 25(t) (7s, 7) • 


In particular, if (16) is satisfied for a curve 70 , then the curve diffusion flow evolving from that 
curve evolves purely by translation. We call the solution a rotator. 

Rotators have been conjectured to exist for the curve shortening flow for some time, with a 
specific conjecture in Altschuler |A . A rigorous classification of these has only recently appeared 
in the literature however Ha . In this section we prove some classification results for rotators. 
Our method is to multiply (16 1 by the curvature and integrate. 

Proposition 14 . Let 7 : S — > R 2 be a smooth closed rotator. Then 7 is a standard round circle. 
Proof: The unit normal vector satisfies 

0 = ^M 2 = {v,v s ) . 

Using the Frenet-Serret frame equations, 

[ n^di = — f k k ss dl = —2 S j (k'ys, 7 ) d£ = 25 f (v s , 7 ) dJL = —2 S f (u, 7 S ) dJL = 0, 

J ■y J -*y J <y J J 

irrespective of S. The curvature of 7 is therefore constant and its image must be a round circle. □ 
Let us now turn to open curves. 

Proposition 15 . Suppose 7 : R —> ]R 2 is a smooth open rotator satisfying 


(17) 


lim 

i —>00 . 


(K(Pi)K s (pi) - 2 S {v{pi),T(pi))) - (K(?i)/c s (?i) - 25 (u(gi),7(gi)) 


= 0 


for arbitrary sequences {pi} and {<7;} asymptotic to —00 and +00 respectively. Then 7(H) is a 
straight line. 
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Proof: Integrate (16) against n to find 

rn 


idl = — 


,dl- 


'Pi 


|9i 
I Pi 


= — 2 S / ( 1 / 5 , 7 ) dl + kk. 


— 2 S (u, 7 ) + KK S 


J Pi 


Taking i —)■ 00 yields the result. 


□ 


Remark: The result of Proposition 15 is false without condition ( [l7| ), for the same reason that 
Proposition 11 is false without condition (12). 


As with translators earlier, we may convert (17) from a pointwise decay condition on the 
curvature to a natural growth condition on n in L z and S ( 177 ) in L [ . 


Proposition 16 . Suppose 7 : 


is a smooth open proper rotator satisfying 


(18) 


lim ( — 


K 2 dl+ 1 -^ 


{v, 7 ) I di = 0 


p-)-oo \p i 7 -i (Sp (0)) p d 7 -i(B 2p (0)) 

where B p ( 0) denotes the closed ball of radius p centred at the origin. Then 7 (R) is a straight 
line. 


Proof: Let r] p = fj p o 7 be the cutoff function used in the proof of Proposition 12 We compute 


K 2 p 4 dl = — 


/ 7 - 1 (s P ( 0)) 


= 2 S 
= -85 


f 

V l (Bp(0) 


KK SS dl — 4 


/ 7 -!(B p ( 0)) 

17 , 7 ) r/ 4 dl — 4 / nn s p s p 3 dl 

l 7 -i(s p (0)) 

u, 7 ) r? s ? 7 3 - 4 / nK s p s p 3 dl 


kk s r] s rf dl 


< £ 


/ 7 - 1 (B p (0)) 

k 2 ?7 4 d£ + 8151 


/ 7 -i(B p (0)) 
4 


/ 7 -i(B p (0)) 


7- 1 (B p (0)) 

*7 7) I M 


k 2 7/ 2 r? 2 df, 


for any e > 0. Choosing £ = 7 and absorbing the first term on the left, we find 

c\S\ 


'dl- 


/ n 2 s p 4 dl< -A- / 

d 7 - 1 (B p ( 0)) P d 7 -i(B 2p (0)) P d 7 - 1 (B 2p (0)) 

Taking p —> 00 in this estimate yields k s = 0 and hence the result. 


u, 7 ) d£. 


□ 
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